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We reply to the comment by Andreotti and Claudin (submitted to Phys. Rev. E, arXiv:0705.3525 ) 
on our paper "Minimal Size of a Barchan Dune" [Phys. Rev. E 75, 011301 (2007)]. We show that 
the equations of the dune model used in our calculations are self-consistent and effectively lead to a 
dependence of the minimal dune size on the wind speed through the saturation length. Furthermore, 
we show that Meridiani Planum ripples are probably not a good reference to estimate the grain size 
of Martian dune sands: the soil in the ripple troughs at the landing site is covered with nonerodible 
elements ("blueberries"), which increase the minimal threshold for saltation by a factor of 2.0. We 
conclude that, in the absence of large fragments as the ones found at the landing site, basaltic grains 
of diameter d = 500 ± 100/im that compose the large, typical dark Martian dunes [K. S. Edgett 
and P. R. Christensen, J. Geophys. Res. 96, 22765 (1991)] probably saltate during the strongest 
storms on Mars. We also show that the wind friction speed u* « 3.0 m/s that we found from 
the calculations of Martian dunes is within the values of maximum wind speeds that occur during 
Martian storms a few times a decade R. E. Arvidson et al, Science 222, 463 (1983); H. J. Moore, 
J. Geophys. Res. 90, 163 (1985); R. Sullivan et al, Nature 436, 58 (2005); D. J. Jerolmack et 
al, J. Geophys. Res. Ill, E12S02 (2006)]. In this manner, the dune model predicts that Martian 
dunes can be formed under present Martian conditions, with no need to assume other conditions of 
wind and atmosphere that could have prevailed in the past. 



In the preceeding comment [l| , Andreotti and Claudin 
pretend to find inconsistencies in the dune model which 
has been used by Parteli et al. [2| in the study of the min- 
imal size of barchan dunes. This model, which consists 
of a coupled set of equations for the wind profile over 
the topography, the sand flux and the evolution of the 
topography with time, has been originally presented in 
refs. later improved in refs. [EH], and repeatedly 

tested through successful quantitative comparison with 
real wind tunnel data and with real dunes measured in 
the field @, [ESI1- In ref. 0, Parteli et al. studied, 
with the dune model, the role of the wind strength and 
inter-dune flux for the shape and the size of the minimal 
dune, and used the results to obtain the wind velocity on 
Mars from the minimal size of Martian dunes. 

The first criticism of Andreotti and Claudin [l[ is that 
the dune model is not self-consistent. They state that the 
saturation length of the sand flux, which determines the 
minimal dune size, should not decrease with the wind ve- 
locity because the relaxation rate is limited by the grain 
inertia. Next, Andreotti and Claudin [l| find that the 
grain size of the ripples at Meridiani Planum landing site 
on Mars is d — 87 ± 25/xm, which is much smaller than 
the grain size d — 500 ± 100 /im of the larger, dark Mar- 
tian dunes, as obtained from thermal inertia data 0] and 
used in the calculations of Parteli et al. [2j. Andreotti 
and Claudin, then, propose an alternative explanation for 
the dependence of the minimal size on the wind speed: 
the effect of slopes. 

The comments of Andreotti and Claudin [j] are con- 
structive and the issues addressed by these authors de- 
serve to be discussed in depth. We organize the present 
reply paper following the same structure of the preced- 



ing comment [l|: Section I, regarding the modeling of 
the flux saturation length and the self-consistency of the 
dune model; Section II, concerning the grain size of Mar- 
tian dunes sand and the reliability of the value of Martian 
wind velocity obtained by Parteli et al. 0; and Section 
III, concerning the effect of slopes on the minimal dune 
size. 



I. SAND TRANSPORT MODEL 

The first criticism of Andreotti and Claudin [l[ refers 
to an apparent inconsistency in our sand transport 
model. They say that, since the grain inertia is not 
included in the evolution of the sand flux, the satura- 
tion length determined by the ejection process can be 
smaller than the length needed for the grains to reach 
their asymptotic trajectory. 

Indeed, in the current model for sand transport we as- 
sume that the characteristic length for the relaxation of 
the mean grain velocity in the saltation layer is much 
smaller than the flux relaxation length determined by 
grain ejection. This can lead to a discrepancy with the 
full model for wind shear velocities it, far from the thresh- 
old Uth- In the following we calculate a modified satu- 
ration length that takes into account both processes and 
show that the saturation length l s is determined by the 
ejection process for the typical range of shear velocities 
found on Earth, i.e. < 3u t h ~ 0.7m/ s. Notice that all 
previous sand dunes simulation results performed with 
the current sand transport model are included within 
both ranges 0, IE 13, IE - 

Following the original approach of Sauermann et al. Q 
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the saltation belt is modeled as a granular fluid layer 
characterized by a vertically averaged mean velocity u 
and grain density p. Both magnitudes obey the mass 
and momentum conservation equations averaged over the 
z-axis. 
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Mass conservation 



The mass conservation over the saltation layer reads 
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where the right hand term accounts for the interchange 
of particles between the saltation layer and the surface 
mainly due to the ejected grains by the splash. This 
term describes the relaxation toward saturation p s of the 
grain density p in the saltation layer. Here, the saturation 
density, defined as the maximum amount of grains carried 
by the wind with a given shear velocity, is given by 
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here U* = u^/uth is the relative wind shear velocity, r t h = 
Pf u th i s threshold shear stress, g is the gravity, 7 and 
a are model parameters and pf is the fluid density [HQ. 



B. Momentum conservation 

Furthermore, the model assumes that the saltation 
layer over a flat surface is only subjected to a mean wind 
drag force and a friction force. The later accounts for the 
momentum lost during the inelastic grain collisions with 
the bed. 

The momentum conservation over a flat bed is given 



by 1 



du 
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where Uf a n is the grain settling velocity. The first right- 
hand term represents a Newtonian drag force exerted by 
the wind with an effective velocity v, while the second 
gives the bed friction. 



II. LINEAR ANALYSIS 

For stationary ID profiles, both conservation equations 
(PandQi reduce to: 
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This coupled system has the homogeneous solution 
(p s ,u s ), that correspond to the saturated state. Intro- 
ducing the linear perturbations around the homogeneous 
solution: p{x) = p s (l + p(x)) and u(x) — u s {\ + u(x)), 
and selecting v(p s ) in the momentum equation, the lin- 
earized system becomes 
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where the characteristic relaxation lengths for the mean 
density and velocity of the saltation layer are respec- 
tively: 
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Here, the mean saturated grain velocity u s is given by 

u s =v s — Uf/V2a (9) 
and the effective wind velocity in the saturated state is 



(TT \ Uth h Zl , Zl ITT 

v s (C/*) = I In 1 (U* 

K \ Zq Z m 



1 



(10) 



where, for simplicity, we select the effective wind velocity 
at the saturated density (v = v s (p s )), and thus the air 
borne shear stress is reduced to the threshold value for 
sand transport 

The relaxation length accounts for the relaxation 
due to the ejection process, while l v includes the grain 
inertia, given by the falling velocity itf a n, in the momen- 
tum balance. 

From the first order linear system Eq. [6] the largest 
relaxation length toward saturation, defined as the satu- 
ration length l s , is given by 
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Figure Q] shows the saturation length l s for Earth con- 
ditions. As was pointed out by Andreotti and Claudin 
the decrease of the saturation length for large shear 
ratios U* is clearly limited by the grain inertia. Indeed, 
we find that the saturation length is determined by the 
ejection process in the range U* < 3.3, which includes 
most real wind conditions. Contrary to our previous as- 
sumptions [1, H, H, @] and in agreement with the sugges- 
tion of Andreotti and Claudin |l| , the spatial relaxation 
of the velocity of the granular layer cannot be neglected 
for larger shear ratios. 
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FIG. 1: Dependence of the density, velocity and overall relax- 
ation lengths, Id, l v and l 3 respectively, with the wind shear 
ratio (Eas, 171 151 and ITTjl for Earth conditions. 



III. THE SIZE AND DENSITY OF GRAINS ON 
MARS 

In the comment, Andreotti and Claudin [jj propose 
that the grains that constitute the sand of Martian dunes 
have diameter d = 87 ± 25 /im. This value has been ob- 
tained by the same authors in a previous work [1 51 ] . in 
which they analysed recent photographs of Martian rip- 
ples taken by the rovers at Meridiani Planum. However, 
the value of grain diameter obtained by Andreotti and 
Claudin [l[ from the analysis of the Meridiani Planum 
ripples is much smaller than the grain size of the typical 
large intra-crater dunes as obtained from thermal iner- 
tia data, i.e. d = 500 ± 100/im [9j. In fact, this value of 
grain diameter, which is in the range of medium to coarse 
sand, has been found, later, to be fairly consistent with 
measurements peformed by several authors in different 
locations on Mars [lj| [TtJ ■ Andreotti and Claudin [l[ , 
thus, conclude that the value of grain diameter obtained 
from the method of thermal inertia, which has been sub- 
ject of research for almost 4 decades, is wrong. 

There is no doubt that the work of Claudin and An- 
dreotti [l5| is of relevance since the measurements of 
grain sizes performed by these authors are based on im- 
ages of unprecedent resolution. However, care must be 
taken before generalizing their results of grain sizes ob- 
tained from the Meridiani Planum ripples to the typical 
large dark dunes on Mars. 



A. Threshold for saltation at Meridiani Planum 

The soil of the Meridiani Planum landing site is covered 
with hematite spherules and fragments reaching milime- 
ters in size. These hematite particles or "blueberries" are 
much larger and denser than the typical basaltic sand 



of Martian dunes. The landing site is in fact a field 
of coarse-grained ripples, whose interiors consist of fine 
basaltic sand in the range of 50 — 125 /im, but which 
are armoured with coarse grains at their crests [HI, . 
Hematite particles with a median diameter of about 1.0 
mm cover more than 75% of the crest area of all ripples. 
On the other hand, the coarse-grain coverage in ripple 
troughs is of almost 50%, whereas inter-ripple areas are 
composed mostly of intact spherules having diameter of 
several milimctcrs, with median 3.0 mm. In comparison, 
"granule" ripple troughs on Earth have an insignificant 
coverage of large fragments, the coarse particles remain- 
ing almost entirely on the ripple crests [19L [20L [2lT | . 

The role of the blueberries for the transport of sand 
at Meridiani Planum ripple troughs is dramatic. It is 
well known that the presence of large particulates shield- 
ing a sand bed increases the minimal wind velocity u^p 
for entrainment of the finer grains into saltation [22j, |23j . 
Gillette and Stockton [22j found experimentally that the 
minimal threshold u*ft of erodible grains with diame- 
ter d ~ 107 /Ltm increased by a factor of k 2.5 in 
the presence of nonerodible grains with diameters about 
D = 2.0 — 4.0 mm having spatial coverage of 45%. In 
fact, the geometrical properties of these experiments are 
very similar to the ones of Meridiani Planum troughs. 
Indeed, a value of k « 2.0 was found later by Nickling 
and McKenna Neuman 23] from experiments with larger 
particles, where d = 270/im and D = 18 mm. 

On the basis of the results from the experiments men- 
tioned in the last paragraph, it was possible to explain the 
formation of the Meridiani Planum ripples. As demon- 
strated in recent publications [H, [lj| , there are strong 
evidences that the minimal wind velocity required to mo- 
bilize the sand grains at Meridiani Planum ripple troughs 
has been effectively increased by a factor k of about 
2.0 — 2.5, as observed in experiments with sand bed 
shielded by nonerodible roughness mentioned above. 

In the absence of nonerodible large fragments, the min- 
imal wind velocity required to entrain sand grai ns into 
saltation can be calculated with the equation [24j |: 



u*ft = A . 



PfLuid)gd 



Pfluid 



(12) 



where g is gravity and pfl u id and p gra i n are the densities of 
the air, respectively of the grains. The Shields parameter 
A is given by (24[: 



A = 0.129 



(1 + 6.0 x 10^ 7 /p grain5 rf 
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,0.5 



A = 0.129(1 + 6.0 X 10" 7 /Pgrain.gd 2 - 5 ) 

■{1 - 0.0858 exp [-0.0617(Re* ft - 10)]} (14) 

for Re*ft > 10, where v is the kinematic viscosity, Re*ft 
is the friction Reynolds number Re*ft = u^d/v, and the 
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constant 6.0 x 10~ 7 has units of kg-m°' 5 -s -2 , while all 
other numbers are dimcnsionlcss. The solid curve in fig. 
[2] shows u*f t as function of the grain diameter calculated 
with Eq. dnj) using p grain = 3200 kg/m 3 , p fluld = 0.02 
kg/m 3 , and v = 6.35 • 10~ 4 m 2 /s 0. 
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FIG. 2: Threshold wind shear velocity ?i*f t for direct particle 
entrainment into saltation on Mars. The solid curve shows 
w*ft calculated with Eq. (|12p . At the Meridiani Planum land- 
ing site, the presence of the blueberries increases it*ft of Mar- 
tian sand by a factor of 2.25 (dashed line). The straight line 
shows the maximum peak of shear velocity, u, = 3.5 m/s, 
associated with the largest storm at the landing site. During 
such a storm, only particles between 69 and 168 /im are ex- 
pected to be entrained by wind at the landing site (dashed 
area). The threshold for the grain diameter of the dark Mar- 
tian sand dunes, d = 500 /im, calculated with Eq. \12), is 
indicated by the filled circle. The empty circle shows u*f t 
for d = 500 /im at the Meridiani Planum landing site: such 
coarse grains cannot saltate at the landing site, at present 
conditions. 



We follow the idea of Jerolmack et al. 19] and cal- 
culate the modified threshold for saltation at Meridiani 
Planum, ku*ft, taking the average value k = 2.25. The 
result is shown by the dashed curve in fig. [21 In this 
figure, the full, straight line represents the maximum al- 
lowed wind friction speed during the gusts of dust storm 
at Meridiani Planum: u* = 3.5 m/s. This value of 
wind speed, which is probably achieved once in inter- 
vals of years [H[ , is estimated to be an upper bound be- 
cause larger wind speeds would result in saltation of the 
hematite spherules, which evidently did not occur dur- 
ing formation of the ripples. As explained previously, the 
winds that formed the ripples at Meridiani Planum land- 
ing site have friction speed u* in the range 2.5 — 3.5 m/s, 
the lower bound corresponding to the minimal threshold 
for creeping motion of the hematite grains [19j . 

Although the estimation of the modified threshold for 
saltation (dashed line of fig. [2]) is very crude, it suggests 
that the wind strength that sculped the soils of Meridiani 
Planum was just sufficient to entrain the grains of small- 
est saltation threshold values, as recognized in ref. [l5l |. 



The dashed area of fig. [2] corresponds to the range of 
grain sizes that are entrained by the wind into saltation 
at Meridiani Planum, 69 < d < 168 /urn, assuming it*f t 
is about 2.25 times the value calculated with Eq. (|12[) . 
The minimum for saltation occurs in fact at about 100 
/im, which is well within the range of grain sizes of the 
sand found in the interior of coarse-grained ripples; on 
the matrix bed in the ripple troughs, and within small 
pits and craters at Meridiani Planum, which apparently 
serve as particle traps [H[ . 

For illustration, the value of u*f t obtained with Eq. 
(fT2"]) for the grain diameter of Martian dunes, d = 500/im, 
is shown by the filled circle in fig. [5] The empty circle 
shows the modified threshold, 2.25u*ft, for d = 500/im. 
We see that saltation of such coarse grains at Meridiani 
Planum landing site would require a wind of u* 5.0 
m/s, which is much larger than the maximum value, 3.5 
m/s. However, it is clear from fig. [2] that basaltic grains 
much larger than those of the landing site can be en- 
trained by a wind of strength u* = 3.5 m/s, in places 
where sand is not shielded by nonerodiblc elements. 

In conclusion, the threshold for saltation transport at 
Meridiani Planum landing site is modified due to the 
presence of nonerodible hematite fragments on the soil. 
Thus, provided other factors as sand induration [25| are 
not affecting the local threshold for saltation, it is very 
plausible that the grains of Martian dunes, which have 
diameter d = 500 ± 100/im [?|, are effectively entrained 
by formative winds of strength 2.5 < u* < 3.5 m/s un- 
der present Martian conditions, since the threshold for 
entrainment of such coarse grains is exceeded at such 
values of u*. 



B. The wind velocity that forms dunes on Mars 

The main criticism in the comment by Andreotti and 
Claudin regarding the results on Martian dunes is that 
the dune model predicts that "very strong" winds [l[ are 
required to form the Martian dunes. We recall the value 
of Martian wind shear velocity obtained in Parteli et al. 
from the minimal dune size [2J: it* « 3.0 m/s. How- 
ever, values of it* about 3.0 m/s are within maximum 
values of shear velocity on Mars, and occur only during 
the strongest dust storms [26|. Sand transport on Mars 
is, thus, expected to consist of short duration events (a 
few minutes) a few times a decade [HI, H3] , and does not 
occur under typical Martian wind velocities, which are 
between 0.3 and 0.7 m/s (28[. We conclude that the value 
it* ~ 3.0 m/s found in Parteli et al. Q from the shape 
of Martian sand dunes is consistent with real values of 
wind velocities expected to occur during sand transport 
on Mars. 

Furthermore, as it will be shown in a future publica- 
tion [3(| , Martian dunes of different shapes and sizes and 
at different locations on Mars can be explained without 
necessity to assume that they were formed "in the past 
under very strong winds" as stated in ref. The cal- 
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culations using the model presented in Parteli et al. @ 
show that the wind velocity on Mars in fact does not 
exceed 3.0 m/s. 

Indeed, the differences in minimal dune sizes on Mars 
can be explained on the basis of different local conditions: 
the average density of the Martian air may vary by a fac- 
tor of 2.0 depending on the location, and in this manner, 
in places where the air is denser (e.g. on the north pole), 
the minimal dune width can be correspondingly smaller 
under the same wind u* — 3.0 m/s. This shear velocity is 
between 1.0 and 2.0 times the threshold for saltation on 
Mars, similarly to the situation on Earth, and is of the 
same order of magnitude as real Martian sand-moving 
winds [lllllllllm 



IV. LINEAR STABILITY ANALYSIS: 
RELATION BETWEEN THE UNSTABLE DUNE 
WAVELENGTH AND THE SATURATION 
LENGTH 

In their comment on "Minimal size of barchan dunes" , 
Andreotti and Claudin 1] also proposed a novel mecha- 
nism to understand the apparent scaling of the minimal 
dune size with the inverse of the wind shear stress, be- 
sides the scaling of the saturation length which arises 
from the derivation of the dune model (Section I). 

Following the work of Rasmussen et al. 31], An- 
dreotti and Claudin include the dependence of the thresh- 
old shear stress u 2 h on the local slope tana = d x h, 
given originally by u 2 h (cosa + sin a/ tan 0) [3l|, where 
9 is the angle of repose, into the linear stability analy- 
sis of the equations for the dune evolution. However, for 
small slopes they approximate the above expression by 
u 2 h (l + d x h/ tan#) instead of the more appropriate one: 
u 2 h (l + (1 — tan8)d x h/ tan#). Furthermore, the scaling 
of the saturated flux with the local slope is not only re- 
duced to the threshold shear stress. In fact, it is more 
complex and can be written as (3lj 



q s = X {ul - u 2 th M)u th /VM (15) 
where M = l + d x h/a, the constant a — tan#/(l — tan (9) 



and u* is the wind shear velocity. 

With this slope dependence the linear perturbation of 
the normalized saturated flux in the Fourier space can be 
written as: 



1 



m-i 



U 2 (A + iB) - — (U 2 + 1) I kh, (16) 
2a ' 1 



where [/» = u*/uth and h is the surface Fourier trans- 
form. 

The linear stability analysis shows that the surface per- 



turbation h is unstable [14[ and has a growth rate a: 



a= i + k'P a (y.) { B -^uf- AMs J (17) 

and thus, the marginally stable wavelength A c , at which 
a = 0, is: 



c ~ U*{2aB-l)-l s[ *' 



(18) 



For a repose angle 6 34° (a w 2) and the typical values 
B = 1.5 and A = 5 [l[, A c becomes: 



A, = S^k(E/*) 
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1 



8tt4 (£/*), 



(19) 



where, due to the factor 2aB — 1 ^> 1, it is clear that the 
main dependence on the shear stress is in the saturation 
length. Notice that if a were simply tan (9 as suggested 
in [l| the constant 2aB 

Summarizing, in the selection of the dominant dune 
wavelength, we find that the "second order" effects in the 
saturated flux apparently counteracts each other. There- 
fore, the contribution of wind speed to the dune size se- 
lection is mainly through the scaling of the saturation 
length, with no simple alternatives. 
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